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ABSTRACT 

We present a study of the morphology, kinematics and ionization structure of the 
extended emission-line regions in four intermediate redshift (0.118 < z < 0.181) ultra- 
luminous infrared galaxies, derived from ARGUS two-dimensional fibre spectroscopy. 

The gas kinematics in the hyperluminous system IRAS F20460+1925 lack coher- 
ent structure, with a FWHM > 1000 kms -1 within 1 arcsec of the nucleus, suggesting 
that any merger is well-advanced. Emission-line intensity ratios point to AGN pho- 
toionization for the excitation of this gas at the systemic velocity. An isolated blob 
~ 8 kpc from the nucleus with a much smaller velocity dispersion may lie in a struc- 
ture similar to the photoionization cones seen in lower-luminosity objects. A second, 
spatially-unresolved, narrow line component is also present on nucleus, blueshifted by 
~ 990 km s -1 from the systemic and plausibly powered by photoionizing shocks. 

IRAS F23060+0505 has more ordered kinematics, with a region of increased 
FWHM coincident with the blue half of a dipolar velocity field. The systemic ve- 
locity rotation curve is asymmetric in appearance, due either to the on-going merger 
or to nuclear dust-obscuration. From a higher-resolution ISIS spectrum, we attribute 
the blue asymmetry in the narrow line profiles to a spatially-resolved nuclear outflow. 
Emission-line intensity ratios suggest shock+precursor ionization for the systemic com- 
ponent, consistent with the X-ray view of a heavily obscured AGN. 

The lower luminosity objects IRAS F01217+0122 and F01003-2238 complete the 
sample. The former has a featureless velocity field with a high FWHM, a high- 
ionization AGN spectrum and a ~ 1 Gyr old starburst continuum. IRAS F01003-2238 
has a dipolar velocity field and an HII-region emission line spectrum with a strong 
blue continuum. After correction for intrinsic extinction, the latter can be reproduced 
with ~ 10 7 05 stars, sufficient to power the bolometric luminosity of the entire galaxy. 

We accommodate this diversity within the merger-induced evolutionary scenario 
for ultraluminous infrared galaxies: the merger status is assessed from the kinematics 
in a way which is consistent with morphological and colour information on the galaxies, 
or with the inferred ages of the young stellar populations and the dominance of the 
AGN. 

Key words: galaxies: individual F01003-2238 - galaxies individual F01217+0122 
galaxies:individual F20460+1925 - galaxies: individual F23060+0505 in- 
frared:galaxies - X-rays: galaxies 



1 INTRODUCTION 

One of the most important results to stem from the IRAS 
mission of the early 1980s was the discovery of a popula- 
tion of galaxies which emit most of their bolometric lumi- 
nosity at mid- to far-infrared wavelengths (see Sanders & 
Mirabel 1996 for a review). Whilst there is little doubt that 
this is thermal re-radiation from heated dust grains, the na- 
ture of the underlying power source remains uncertain - for 
some individual objects and for the class of ultraluminous 



infrared galaxies (ULIRGs, with Lib. > 10 12 Lq) as a whole. 
The resolution of this problem has significant implications 
for the wider context of extragalactic astronomy. For exam- 
ple, some of the most luminous ULIRGs and some hyper- 
luminous objects (I/ir > 10 13 Lq) exhibit high-ionization 
Seyfert-2-like optical spectra [eg. IRAS F15307+3252, Hines 
et al. (1995); IRAS P09104+4109, Hines & Wills (1993)], 
leading to suggestions that they constitute the long-sought 
'type-2' quasars - active galactic nuclei (AGN) with Seyfert 
2-like optical spectra but at QSO luminosities. The results 
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of mid-infrared ISO spectroscopic surveys (Genzel et al. and 
Lutz et al. 1998), however, suggest that more than 80 per 
cent of ULIRGs are predominantly powered by recently- 
formed massive stars, with the AGN-powered fraction in- 
creasing to around half above an infrared luminosity of 
2 x 10 12 L Q . 

For the past decade, the formation and evolution of 
ULIRGs have been discussed in the context of the paradigm 
of Sanders et al. (1988a), in which they constitute an early 
phase in the merger-induced formation of optically-selected 
quasars. Following the merger of two gas- rich spiral galaxies, 
tidal torques drive material into the merger nuclei, leading 
to prodigious star formation and the fuelling of supermas- 
sive black holes. The circumnuclear regions are subsequently 
swept clear of obscuring dust by radiation pressure, exposing 
the quasar nucleus. Whilst this model was in part suggested 
by the frequency of occurrence of ULIRGs in interacting or 
merging systems (Sanders et al. 1988b), since corroborated 
in the study of larger samples (eg. Clements et al. 1996) , it 
is undoubtedly not the whole story. Indeed, from their ISO 
spectroscopy, Genzel et al. (1998) found no obvious trend 
for the AGN component to dominate in the more advanced 
mergers of their sample, implying that merger-phase alone 
does not determine the relative dominance of an AGN com- 
ponent. Similarly, the Fabry-Perot Ha velocity mapping of 
four ULIRGs by Mihos & Bothun (1998) showed that ul- 
traluminous activity is not confined to late stage mergers, 
suggesting that details such as the internal structure or gas 
content of the merging galaxies can strongly influence their 
evolution. 

In a study of the extended emission line region (EELR) 
of the ULIRG IRAS P09104+4109 with the ARGUS integral 
field spectrograph on the Canada France Hawaii Telescope 
(CFHT), Crawford & Vanderriest (1996) found that the gas 
was energised by an obscured quasar continuum and sug- 
gested that a 1200 kms -1 nuclear outflow was responsible 
for displacing X-ray-emitting gas from the centre of the sur- 
rounding cooling-flow cluster. These findings prompted us to 
use ARGUS to study a small sample of ULIRGs, in order to 
discover whether the properties of IRAS P09104+4109 are 
typical of such luminous objects or due to its unusual envi- 
ronment, and whether any trends with increasing luminosity 
can be identified. 

The two principal targets in our new sample are 
the famous QSO-like systems IRAS F20460+1925 and 
F23060+0505, at redshifts of 0.181 and 0.174 respec- 
tively; they were chosen partly for their large luminosities 
(falling just either side of the conventionally-adopted ul- 
tra/hyperluminous divide) and also because they have been 
studied extensively at other wavelengths, with evidence from 
both X-ray studies and near-IR broad emission lines for ob- 
scured active nuclei. Included for comparison are two lower- 
luminosity, lower-redshift objects, IRAS F01217+0122 and 
F01003-2238, the latter having Wolf-Rayet features in its 
spectrum, suggestive of a decaying starburst (Armus, Heck- 
man & Miley 1988). The primary aim of the present area- 
spectroscopic study is to understand the observed kinemat- 
ics and morphology of these EELRs within the context of 
plausible dynamical models for their origin in, for exam- 
ple, galaxy mergers. Secondly, we seek to assess the viability 
of various ionization mechanisms, including power-law pho- 



toionization, fast shocks and their precursor HII regions, and 
young massive stars. 

For each of the two most luminous objects, we discuss 
the views provided by the ARGUS data on the gas kine- 
matics and its ionization structure. We also present higher 
resolution long-slit spectra of these two objects taken with 
the William Herschel Telescope (WHT), and comment upon 
multicolour images from the Jacobus-Kapteyn Telescope 
(JKT). In obtaining the former, we sought both to confirm 
the ARGUS findings and to overcome some of their inherent 
limitations arising from their modest spectral resolution and 
decreased signal-to-noise towards the blue part of the spec- 
trum. Following a less detailed presentation of the ARGUS 
data for the other two objects, we conclude with a discus- 
sion of the comparative properties of the sample as a whole. 
For convenience, we summarise in Table 1 the lay-out of the 
results section. 

The cosmological parameters Ho = 50kms _1 Mpc -1 
and qo = 0.5 have been adopted throughout. 

2 OBSERVATIONS, DATA REDUCTION & 
ANALYSIS 

The details of all the observations are presented in Table 2 
and discussed separately below. 

2.1 ARGUS integral field spectroscopy 

Observations were made with the ARGUS integral field 
spectrograph on the CFHT during the nights of 1996 August 
14 and 15. The instrument consists of a 12.8 x 7.8 arcsec 2 
aperture in the focal plane of the telescope, within which 
0.4 arcsec diameter optical fibres are arranged into 27 rows. 
At the entrance to the multi-object spectrograph the fibre 
bundle is drawn out to mimic a long slit arrangement with 
adjacent rows being separated by blank fibres. The O300 
grism and STIS 2 CCD were employed, covering the range 
4000-10000A at a dispersion of 5A per pixel and a resolution 
of approximately 10. 3A FWHM. A more complete descrip- 
tion of the instrumental set-up and performance is provided 
by Crawford and Vanderriest (1997). 

The principal targets, IRAS F20460+1925 and 
F23060+0505, were observed for total exposures of 10800 
and 9000 seconds, respectively. The seeing during this pe- 
riod was typically 0.5-0.7 arcsec FWHM, but it deteriorated 
three hours into the second night to ~ 1.5 — 2 arcsec FWHM 
when the lower luminosity objects, IRAS F01003-2238 and 
F01217+0122, were observed (for 3000 and 5400 seconds, 
respectively). 

The reduction of the ARGUS data was performed 
within IRAF. Cosmic rays were removed by co-adding the 
separate images taken of each field. Thereafter, a two- 
dimensional transformation was applied to the CCD frames 
in order to correct for a form of optical distortion in which 
the blank spacer fibres were seen to bow outwards towards 
the left- and right-hand edges of the frame. Each fibre spec- 
trum then projected onto two rows of the chip across the 
wavelength range. Wavelength calibration was performed by 
applying a wavelength-pixel number transformation derived 
from He/Ne, He/Ar, He/Ne/Ar and Ar arc spectra. Bias 
subtraction was problematic, as the level varied throughout 
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Table 1. Lay-out of the results sections of this paper 
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Table 2. Summary of the observations 
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ARGUS integral field spectroscopy on the CFHT 






1996 August 14-15 








Spectral resolution: 10.3 A FWHM 






Wavelength coverage: 4000-10000 A 




Seeing: arcscc 
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Integration time: s 
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ISIS longslit spectroscopy on the WHT 






1998 September 21 








Spectral resolution: 3.1 A FWHM 






Wavelength coverage: 3650-6000 
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Seeing: (arcsec) 
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Integration time: (s) 




3600 
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JKT multi-band imaging 








1998 August 20-21 






Seeing: (arcsec) 




~ 1 


~ 1 


B-band Integration time: (s) 




2000 


2000 


R-band Integration time: (s) 




1200 


1200 


I-band Integration time: (s) 




2000 


1800 



t ^FIR = £(40 — 500/irra), computed from IRAS fluxes using the formula in Sanders & Mirabel (1996). 

\ -L Bol computed by Frogel et al. (1989) 
* L(0.5 - 120/oti), given by Hill et al. (1991) 



each night and with spatial location on the chip. In conse- 
quence, a separate bias frame was synthesised for each data 
frame as follows: all of the good bias frames with similar 
(x,y) structure were combined to produce a normalised bias 
frame, which was then multiplied by the y-variation of the 
x-averaged bias level for the data frame under consideration 
(as estimated from the blue region of the chip where the 
quantum efficiency of the system is low). Flat-fielding was 
not possible due to features in the flat-field exposures (pos- 
sibly induced by local heating of the CCD) which rendered 
impossible their bias subtraction according to the above pro- 
cedure. Compensation for spectrally-averaged fibre-to-fibre 
sensitivity variations was made at a later stage of the anal- 
ysis. This was done by using an arc-lamp exposure to select 
those fibres with a wavelength-integrated flux reduced by 
more than ~ 10 per cent with respect to adjacent fibres. 
The derived transmission correction factors for such fibres 



were later used to scale their observed emission lines fluxes. 
This local method does not correct for global variations in 
the sensitivity of the system across the aperture, as they 
cannot be distinguished from any spatial gradient in the arc 
light illumination. 

Although spectra were taken of six flux standard stars, 
it was subsequently found that only two have fluxes tab- 
ulated up to 10000A (the rest terminate at 8400A) and, 
furthermore, that one of these is contaminated by an un- 
physical feature due to a second order of interference. We 
thus had only one standard star covering the spectral range 
of interest. At wavelengths less than 8400 A a flux uncer- 
tainty of < 25 per cent is estimated, based upon the scatter 
between the flux sensitivity curves in this region. Correc- 
tions were made for atmospheric extinction and Galactic 
reddening, assuming for the latter the A*h : E(B — V) con- 
version of Bohlin, Savage & Drake (1978) and an R value of 



© 0000 RAS, MNRAS 000, 000-000 



4 



3.2 (line-of-sight Galactic hydrogen densities were obtained 
from Stark et al. 1992). The CCD frames were split into 
individual fibre spectra using an ARGUS-dedicated 1RAF 
package. Sky subtraction used an average sky spectrum for 
each row of the hexagon, constructed from the mean of four 
or five fibre spectra in each of that row and the adjacent 
two. Care was taken to avoid fibres in which the effects of 
aperture vignetting were not negligible. 

Comparison of the emission-line fluxes derived from 
the ARGUS data with previously-published values indicates 
that our fluxes are systematically larger, by at least 50 per 
cent. This may be partly due to the fact that previous 
long-slit spectra did not admit all of the extended emis- 
sion. Comparison of the [OIII] and [Oil] fluxes with those 
of Hough ct al. (1991) for IRAS F23060+0505 and with the 
ISIS data (section 2.2) confirms the reliability of the relative 
flux calibration of the ARGUS data blueward of redshifted 
[OIIIJA5007 at ~ 5900A. The ARGUS flux calibration does 
not vary across the aperture, and errors in its absolute level 
do not significantly affect the scientific conclusions of this 
paper. The latter depend mainly upon the kinematic de- 
composition of individual line profiles and upon flux ratios 
of lines either close in wavelength (eg. [OIII]/H/3) or from 
spectral regions over which the relative flux calibration is 
reliable (eg. [OIII] /[Oil]). 

2.2 ISIS spectroscopy 

Higher-resolution long-slit spectra of IRAS F20460+1925 
and F23060+0505 were obtained with the blue arm of the 
ISIS spectrograph on the 4.2m WHT during the night of 
1998 September 21. A one hour integration was performed 
for each object with the R300B grating in place. In con- 
junction with the EEV12 CCD this yielded a wavelength 
coverage of 3650-6500A at a spatial scale of 0.2 arcsec per 
pixel and a spectral resolution of 3.lA FWHM. The width 
of the slit was matched to the seeing (~ 1 arcsec) and ori- 
ented at a position angle chosen to pass through features of 
potential interest (described later). 

The data reduction was performed within IRAF. Fol- 
lowing bias subtraction, the data frame was rotated by 
~ 0.3 degrees to bring the slit and dispersion axes into more 
exact alignment with those of the chip. The spectra were 
subsequently wavelength-calibrated using a CuAr arclamp 
exposure, flux-calibrated using observations of the standard 
star HD192281, corrected for atmospheric extinction and 
Galactic reddening, and sky-subtracted. Since no structure 
could be discerned in the flat field frames, the data were not 
flat-fielded. 



2.3 JKT multicolour imaging 

Using the lm Jacobus-Kapteyn Telescope (JKT), we ob- 
tained images of IRAS F20460+1925 and F23060+0505 in 
the B, R and I bands. The observations were carried out 
on the nights 1998 August 20-21, with integration times of 
between 1200 and 2000 seconds per object, per filter, in see- 
ing conditions of approximately 1 arcsec FWHM. The Tek4 
CCD was employed, covering a 5.6 square arcmin field of 
view on a scale of 0.33 arcsec per pixel. 

The preliminary reduction of these data, namely bias 



and sky-subtraction, was performed within IRAF, but they 
were not considered to be of sufficient quality to merit flux- 
calibration and detailed quantitative analysis. 

2.4 Emission-line fitting: some general remarks 

The most prominent emission line complexes in the ARGUS 
spectra of all five objects are those of [OIII] A4959, 5007+H/3 
and Hq+[NII]A6548,6584. The [OIIJA3727 doublet is de- 
tected only in the nuclear fibres as the signal-to-noise ra- 
tio deteriorates significantly towards the blue end of the 
spectrum. Due to deblending uncertainties in the [NII]+Ha 
complex of these relatively low resolution spectra, kinematic 
information was extracted only from fits to the [OIII]+H/3 
profile. Fitting was performed using the QDP package (Ten- 
nant 1991). 

In the first instance, the narrow lines in a complex were 
fitted with single component gaussian profiles having the 
same velocity width and redshift. The [OIII]A5007/4959 flux 
ratio was set to the theoretical value of 3:1 and the contin- 
uum was modelled as a constant. The fitting of such a profile 
facilitated the rapid identification of the gross kinematical 
features, providing some indication of those aspects worthy 
of more diligent fitting with more complex models. Many of 
the emission lines were found to be asymmetric in appear- 
ance and thus poorly described by a single gaussian. In most 
such cases the total line profiles could be satisfactorily mod- 
elled by a function consisting of two [OIII]+H/3 complexes 
of the above form with a fixed velocity offset between them, 
and with the [OIII] and H/3 fluxes of the second compo- 
nent expressed as fractions of those in the first. The velocity 
widths of all lines in the two components were, however, con- 
strained to be equal. Note that in fitting a double component 
profile we do not wish to imply that there are precisely two 
physically-distinct components to the emission. Asymmetric 
and multi-component line profiles are a common feature of 
the narrow line spectra of Seyfert galaxies and a discussion 
of the complexities associated with their physical interpre- 
tation is given by, for example, Veilleux (1991). 

The decomposition of the Ha+[NII] complex was not 
always possible, because the presence of multiple veloc- 
ity components, and in some cases broad Ha, frequently 
rendered the entire feature unresolved. Where it could 
be modelled with reasonable confidence, multiple velocity 
component profiles were employed, constructed in exact 
analogy to those described above for the [OIII]+H/3 com- 
plex. Specifically, the lines within both velocity components 
were assumed to have identical velocity widths and the 
[NII]A6584/6548 flux ratio was set to the theoretical value of 
3:1. Any broad Ha was assumed to be at the redshift of the 
systemic narrow line component. The [OII]A3727 doublet 
was modelled as a single gaussian, which is a satisfactory 
approximation at this resolution. 

In order to understand how systematic effects might 
affect the interpretation of the parameter values extracted 
from fits with simple models, we performed some simula- 
tions. Synthetic, noise-free, profiles for the [OIII]+H/3 com- 
plex were generated and fitted with the above double compo- 
nent model. The findings are alluded to at various points in 
the text, and discussed in detail in Appendix A. The higher 
resolution of the ISIS spectra allowed us to fit the [OIII]+H/3 
complex with a model in which the velocity widths of the 
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two components were not constrained to be the same, thus 
overcoming these systematic effects. The ISIS data also have 
significantly better wavelength coverage and signal-to-noise 
blueward of H/3 than the ARGUS data, thus permitting 
the analysis of emission lines such as those of [OIIJA3727, 
[Nelll], [NeV], Hell and higher order Balmer lines. Where 
two (as opposed to just a single) gaussian velocity compo- 
nents could justifiably be fitted to such profiles, their red- 
shifts and widths were fixed at the values derived from the 
[OIII]+H/3 complex. 

The position-dependent instrumental resolution was 
subtracted in quadrature from the fitted emission-line 
widths, with the former determined from fits to groups of 
arc lines at wavelengths similar to the galactic emission-lines 
under consideration. All FWHM quoted in this paper have 
thus been corrected for this instrumental effect. 

There are two additional factors which serve to compli- 
cate the profile fitting and its subsequent interpretation. The 
first is that projection effects result in the observed profile 
at any point in the sky being due to emission from all depths 
along the line of sight. The second is a consequence of the fi- 
nite size of the point spread function (PSF) of the telescope 
optics convolved with the atmospheric seeing, as sampled 
by the fibre arrangement: the light from a point source is 
collected by several fibres. As a result, the profile observed 
in a given fibre is contaminated by light from points on the 
sky centred on neighbouring fibres. A proper understanding 
of this effect is necessary when the observed object consists 
of a bright point source enveloped by low surface brightness 
diffuse emission, as is the case for several of the ULIRGs in 
our sample. The standard star exposures taken at frequent 
intervals enable the structure of the PSF to be monitored: 
its FWHM was typically found to be larger along the rows 
of fibres than in the perpendicular direction (see Crawford 
& Vanderriest 1996). 

All errors quoted in this paper are A\ 2 = 2.7 intervals 
for one parameter of interest (corresponding to 90 per cent 
confidence), assuming that the original fit was good. They 
are fitting errors only and do not incorporate any systematic 
effects from the data reduction. 



3 RESULTS FOR IRAS F20460+1925 

IRAS F20460+1925 was selected from the IRAS catalogue 
by Frogel et al. (1989) as part of a programme to study 
sources with flat far-infrared spectral energy distributions. 
They classified it as a Seyfert 2 galaxy at z=0.181 with ex- 
ceptionally broad permitted and forbidden lines (FWHM be- 
tween 1700 and 2100 kms" 1 ). Hines (1991) and Veilleux et 
al. (1997) detected broad Paa with FWHM ~ 3000 km s -1 
and deduced intrinsic extinctions to the hidden BLR of 
Ay > 3.6 and > 6.4 mag, respectively (cf. the Av = 2.4 
mag to the NLR found by Frogel et al.). These values may 
not represent that along the line of sight to the BLR if a 
substantial portion of the BLR spectrum is contributed by 
scattered light, as is suggested by Kay & Miller's (1989) 
detection of broad H/3 in polarised light. 

Ogasaka et al. (1997) reached a similar conclusion from 
their 2-10 keV ASCA spectrum, which they modelled as 
a power-law with photon index F ~ 2 along with iron 
K emission and excess photoelectric absorption of Nh ~ 



2.6 x 10 22 cm -2 . They deduced, however, an absorption- 
corrected value of Lx/£bo1 an order of magnitude below 
that of a typical AGN, perhaps because the direct emis- 
sion is completely blocked by an obscuring torus; that even 
the scattered X-rays are then absorbed indicates that the 
nucleus of this galaxy is indeed heavily obscured. Alterna- 
tively, it could be that a starburst contributes significantly 
to the bolometric luminosity, but their data neither support 
nor refute the existence of an additional soft thermal com- 
ponent. 

Young et al. (1996a) included the object in their sample 
of "warm" IRAS galaxies with Seyfert- 2-like optical spectra 
for which they presented polarimetry and modelling. They 
found multicomponent [OIHJ+H/? lines and an [NII]+Ha 
complex which could be approximated by double component 
narrow lines along with broad Ha. Based on measurements 
of the broad and polarised Ha fluxes, they deduced an in- 
trinsic polarisation of 12 per cent and an extinction of 1.2 
mag in the scattered light. Along the line of sight to the 
near-IR-emitting regions they calculated an Av of 14 mag, 
consistent with the lower limit on the extinction to the BLR 
deduced by Hines (as the BLR is expected to lie within the 
near-IR-emitting region) . 



3.1 The ARGUS data 

3.1.1 Kinematics and spatial distribution of the 
line-emitting gas 

An examination of the CCD frames for this object prior 
to their decomposition into fibre spectra clearly showed the 
double component nature of the [OIII]+H/3 complex in and 
around the nucleus. A localised region of much narrower, 
single component emission, spatially isolated from the nu- 
cleus and redshifted relative to it, could also be discerned 
and will henceforth be referred to as the north-eastern blob. 
Fig. [l] shows the [OIIIJ+H/3 complex of the nuclear fibre, the 
north-eastern blob and an additional kinematically distinct 
area to be referred to later as the western blob. 

When the spectra were fitted with the double compo- 
nent [OIII]+H/3 model, it was found that two components 
were present out to a radius of ~ 2 arcsec from the fi- 
bre of peak intensity (the nucleus) [1 arcsec = 4kpc at 
z=0.181]. In the nucleus itself, the blueshifted component 
is offset from the systemic by ~ — 975 km s -1 with relative 
intensities in [OIII] and H/3 of ~ 1.6 and 0.9 respectively, 
implying that it is the more highly ionized of the two com- 
ponents. Fig. |^ shows reconstructed images of the object 
in several (continuum-subtracted) emission lines. [OII]A3727 
could be detected only within a limited region and reliably 
decomposed into two velocity components over a yet smaller 
area, so the flux in the [Oil] map is that extracted from a 
fit to a single gaussian. A rigorous fitting of the [NII]+Ha 
complex could not be performed because of the difficulty 
of deblending the two (roughly equal-intensity) components 
to each line and the broad Ha; we thus map the flux in a 
(continuum-subtracted) cut between 7730 and 7780 A, wide 
enough to include the broad Ha. The fluxes in all of the lines 
peak on the nuclear fibre (marked by a cross) for both the 
systemic and blueshifted components and exhibit the same, 
essentially circular, morphology. The total luminosity of the 
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Figure 1. The [OIII]+H/3 complex in the nuclear fibre, and 
the kinematically distinct north-eastern and western blobs 
in IRAS F20460+1925. The position of the nuclear fibre sys- 
temic velocity is indicated by the dashed line. 



See separate figure 



Figure 2. Reconstructed images of IRAS F20460+1925 in 
various emission lines, as described in the text. The square 
root of the integrated line flux has been plotted on a com- 
mon scale (in units of 10 _8 erg ' 5 cm _1 s -0 ' 5 ) in order to en- 
hance the dynamic range; the north-eastern blob is high- 
lighted with diamonds. Discs are drawn only for those fibres 
with a significant detection. The actual flux ranges spanned 
by the fibres shown arc (in units of 10 -15 erg cm -2 s _1 ): 
0.021-8.53 (for systemic [OIII]), 0.063-14.0 (for blueshifted 
[OIII]), 0.16-8.26 (for [Oil]), and 0.063-7.48 (for [NII]+H«). 
All maps peak on the same fibre, the nucleus, marked with 
a cross. North is at the top and east to the left in this figure. 

system in the [OIII]A5007 line is 2.56 x 10 43 ergs" 1 , 54 per 
cent of which comes from the blueshifted component. 

Although a blueshifted component to the [OIII] emis- 
sion is detected in fibres up to 2 arcsec from the nucleus, 
this may not necessarily indicate an outflow with a radius 
of 8kpc. Caution is necessary because the morphology of 
the second component superficially resembles that of the 



telescope PSF and thus might in fact be unresolved. To 
investigate this, a PSF for the wavelength interval 5500- 
6000A was constructed from the observation of the flux 
standard BD404032 taken immediately after the first of the 
IRAS F20460+1925 fields (and thus under comparable see- 
ing). The intensity of the blueshifted component (relative 
to that in the nuclear fibre) was measured in eight fibres 
at various distances from the nucleus and compared with 
that expected from a point source. Within the sampling un- 
certainties arising from the fact that the standard star PSF 
appears to peak between fibres whereas all the intensity maps 
in Fig. |^ peak sharply on a single fibre (to say nothing of the 
profile fitting uncertainties) , the second component appears 
unresolved. A further argument in support of this conclu- 
sion is that the line-of-sight velocity of an outflow with any 
reasonable geometry would exhibit spatial structure (eg. de- 
creasing to zero at the edges for a spherical bubble) , whereas 
the observed wavelength of the blue component is constant 
(to within the profile- fitting uncertainties) in all fibres where 
it is seen. 

Fig. ^ maps the velocity and FWHM of the systemic 
[OIII]+H/3 component. An asterisk on a fibre indicates that 
the fitted width is consistent with the profile being spec- 
trally unresolved and for such fibres, the plotted FWHM is 
taken as the average of that in neighbouring fibres. From 
this figure, it can be seen that the velocity of the systemic 
component is substantially constant across the nebula, apart 
from two kinematically distinct regions. The north-eastern 
blob is centred at a projected radius of 8kpc, with single 
component emission redshifted by ~ 240 km s" 1 relative to 
the nucleus, and the western blob, an approximately ellipti- 
cal region to the west of the nucleus, is at a relative velocity 
of ~ 400kms _1 . Although the velocity separation of the 
two components is poorly-constrained in the western blob, 
it is several hundred kms -1 larger than it is in the nucleus. 
As the latter figure is comparable to the difference in the 
systemic velocities of the nucleus and the blob itself, this is 
consistent with the blueshifted component being unresolved. 

It can be seen in Fig. ^ that, with the exception of fi- 
bres in the north-eastern blob wherein just a single narrow 
(and in places unresolved) component of emission is seen, 
the lines have a FWHM of several hundred kms - , exceed- 
ing 1000 km s" 1 within ~ 1 arcsec of the nucleus. Such large 
velocity dispersions suggest that the gas is not in a disk-like 
structure and that it has significant motions in all three spa- 
tial dimensions. This situation is to be contrasted with the 
findings of Unger et al. (1987) who, in their study of the 
extended narrow line regions of Seyfert 2 galaxies, deduced 
that the kinematics of the gas within the ionization cones 
of such objects were consistent with normal galactic rota- 
tion. Typical velocity dispersions for neutral gas in the discs 
of normal spiral galaxies are 7-10 kms -1 (van der Kruit & 
Shostak 1984). For IRAS F20460+1925 it is remarkable that, 
on kiloparsec scales, the FWHM exceeds the upper end of 
the range (300 — 500 kms -1 ) found for classical narrow line 
regions in AGN, which are typically 10 — 1000 pc in size. 

There is a caveat to our interpretation of the fitted 
velocity field. Due to the relatively low resolution of the 
ARGUS spectrograph, the systemic and blueshifted compo- 
nents were constrained to have the same width. It was, how- 
ever, discovered during emission-line fitting simulations (see 
Appendix A) that the parameters extracted from such fits 
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Figure 3. The upper map shows the velocity width FWHM 
common to the systemic and blueshifted components in 
a two-component fit to the [OIII]+H/3 complex of IRAS 
F20460+1925, corrected for the position-dependent instru- 
mental resolution. Those fibres marked with asterisks have 
an unresolved linewidth (see text). The lower map shows 
the radial velocity of the systemic component of [OIIIJ+H/3 
emission relative to the nuclear fibre (marked with a cross) . 
North is at the top and east to the left in this figure, and 
the kinematically distinct north-eastern and western blobs 
are indicated. 

systematically mis-represent the true properties of the two 
components (especially their relative intensities and veloc- 
ity separation) if their widths are in fact markedly different. 
This qualifies our deduction that the second component is 
spatially unresolved, and reduces the physical significance 
which should be attached to the magnitude of its measured 
velocity offset from the systemic. That the blueshifted com- 
ponent appears to be spatially unresolved implies that its 
profile results from integration over an entire velocity field 
and thus suggests that it may be much broader than the 
systemic. The resolution of these potential sources of ambi- 
guity was one of the motivations behind our obtaining the 
higher signal-to-noise ISIS spectrum of this object, which is 
discussed in section 3.2. 

3.1.2 Ionization structure 

The ionization structure of the gas has been derived from 
the [OIII] /H/3 and [OIII]/[OII] ratios, since the analysis of 
the [NII]+Ha complex is plagued by deblending uncertain- 
ties. Whilst [OIII]A5007 is sufficiently strong to permit its 
use as a kinematic diagnostic over the region covered by 
the maps in Fig. ^, the H/3 flux in the individual fibres is 
poorly constrained beyond radii of ~ 1 arcsec and is thus 
measured in groups. Fig. ^ maps the (systemic) [OIII]/H/3 
and [OIII]/[OII] ratios. Due to the low signal-to-noise the 
[OII]A3727 flux used here is taken from a fit with a single 
gaussian and, for consistency, that of [OIII] is the sum of 
the systemic and blueshifted components. All emission-line 
ratios are the observed values (ie. without correction for in- 
trinsic extinction) unless stated otherwise. 

In Fig. ^, the line ratios from these fibres are shown 
along with points representing the blueshifted component 
(as determined from a fit to a group of seven fibres centred 
on the nucleus) , the kinematically distinct north-eastern and 
western blobs, and two other regions of the extended emis- 
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Figure 4. The up- 

per figure maps the ratio of [OIII] A5007/ [Oil] A3727 in all 
fibres of IRAS F20460+1925 where [Oil] is detected, for the 
sum of systemic and blueshifted components (see text). The 
lower figure maps the variation of [OIII] /H/3 in the systemic 
component; H/3 could be detected only in the individual fi- 
bres shown outlined; elsewhere, fits were carried out to the 
spectra of groups of fibres (not outlined). North is to the top 
and east to the left. 




[OII]/[OIII] 

Figure 5. An ionization state diagnostic diagram for IRAS 
F20460+1925. The open circles represent the central seven 
nuclear fibres, filled circles individual fibres around the nu- 
cleus with detectable H/3, the filled star the blueshifted com- 
ponent, the filled box the north-eastern blob and the filled 
triangles other groups of off-nucleus fibres, eg. the western 
blob. The dashed box defines the shock+precursor locus of 
Dopita & Sutherland (1995) and the continuous lines delin- 
eate a series of CLOUDY photoionization models (see text 
for details). The plotted ratios are those observed, ie. with- 
out correction for intrinsic reddening; the arrow denotes the 
correction factor for an intrinsic extinction of Av ~ 2 mag 
(see section 4.2). The dotted cross represents a typical size 
for the error bars on each point, which have been omitted 
for clarity. 
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sion to the south-east and north-west of the nucleus, re- 
spectively. Except for the point representing the blueshifted 
component, the systemic [OIII] and H/3 fluxes are used to 
calculate the line ratios, as are the [Oil] fluxes for the seven 
nuclear fibres, in which decomposition into two velocity com- 
ponents was possible. It can be seen that the blueshifted 
component is more highly ionized than the systemic, and 
that the north-eastern blob is also at a high excitation level. 

The purpose of producing Fig. ^ is to use it in con- 
junction with the results of ionization models to understand 
how the line emission is powered and also to examine how 
the ionization state may vary across the system. The two 
continuous lines on this figure were obtained from a grid of 
CLOUDY (version 90, Ferland 1996) models using the 'table 
agn' continuum. They represent models in which the hydro- 
gen nucleon density at the exposed face of the (ionization- 
bounded) clouds assumes values of 10°' 5 and 10 35 cm , 
respectively, and along which the ionization parameter U 
varies between 10~ 3 ' 5 and 1CT 1 . The 'table agn' continuum 
in CLOUDY is similar to that observed by Mathews and Fer- 
land (1987) for a typical radio quiet active galaxy, and as 
such it may differ from any photoionizing continuum in this 
object, concerning the true shape of which we are largely ig- 
norant. At the higher density, the models stopped due to the 
build-up of radiation pressure for U > 10 -2 ' 75 and are thus 
not shown in Fig. ^| The dashed lines encompass the predic- 
tions of the shock + prescursor ionization models of Dopita 
and Sutherland (1995), with the shock velocity increasing 
from 200 to 500 kms -1 vertically and magnetic parameter 
{B/n - 5 ) in the range 0-4 G cm" 15 . 

Whilst a single diagnostic diagram cannot be used to 
discriminate between the two ionization mechanisms, Fig. ^ 
does not even suggest the dominance of one over the other: it 
could be that there is a mixture of unresolved velocity com- 
ponents at different ionization levels. Indeed, pure photoion- 
ization by an AGN continuum seems unlikely in view of X- 
ray data which show that the AGN continuum source is ob- 
scured by a column density of at least Nh — 2.6 x 10 22 cm -2 
(Ogasaka et al. 1997); an optical- UV continuum of the req- 
uisite strength could not emerge through such obscuration. 
In addition, if the measured line-widths are entirely kine- 
matic in origin, shocks would be expected to be an important 
source of ionization. It should, however, be noted that the 
shock models of Dopita & Sutherland do not extend to veloc- 
ities in excess of 500 kms -1 , whereas we measure a FWHM 
of around twice this value. No obvious correlation between 
the kinematics of the gas and its ionization state (such as 
would exist for shock excitation) is apparent. Neither is it 
possible at this signal-to-noise to measure other shock di- 
agnostic line ratios (eg. HeIIA4686/H/3 or [NeV]A3426/H/3) 
and thus to self-consistently implicate shock excitation. All 
that can be deduced with confidence from the ARGUS data 
is that ordinary O and B stars are not an important source 
of ionization. 

We comment here on the appearance of the [OIII] /[Oil] 
map in Fig. ^. The increase in the ratio north-east of the 
nucleus may result from the escape of nuclear continuum 
radiation along a relatively unobscured line of sight. It is 
interesting to speculate on the association of this structure 
with the north-eastern blob, which does not appear on the 
[OIII] /[Oil] map of Fig. ^ (a 3a lower-limit on this quantity 
therein is 2.5), but which stands out clearly as a region of 



high excitation in [OIHJ/H/3 (= 12 ± 4) at the same posi- 
tion angle. We wish to determine whether the north-eastern 
blob could be powered by continuum radiation from the ac- 
tive nucleus, or whether some in situ ionization source is 
required. 

To address this question, we compare the rate at which 
the nuclear continuum supplies hydrogen-ionizing photons 
to the blob, to the rates required by CLOUDY models 
which are able to account for its observed [OIII]/H/3 value 
and dereddened [OIII]A5007 luminosity (1.34 x 10 42 ergs" 1 ). 
For the nuclear continuum we extrapolate the X-ray spec- 
trum found by Ogasaka et al. (1997) with a power-law of 
photon index 2. We do not normalize it to the value of 
2.7 x 10 44 ergs _1 which the latter authors deduced for the 
absorption-corrected 2-10 keV luminosity, because, due to 
scattering, the true nuclear luminosity may be some un- 
known factor, fscatt, larger. Assuming the blob to lie at its 
projected radius of approximately 8kpc and to present an 
area A c i ou d kpc 2 of exposed cloud-face to the incident con- 
tinuum, we integrate the power-law between Ai = 228A and 
A2 = 912A (ie. between the Lyman limit and the energy 
beyond which photons are used to ionize Hell rather than 
HI), to deduce that the active nucleus supplies the blob 
with hydrogen-ionizing photons at a rate AT.' 1 ^ = 7.2 x 
10 51 Adoudfscatt s -1 . For given values of fscatt and A c i oud , 
the latter estimate of is likely to be lower limit because 
the true spectrum may have a blue-bump which exceeds the 
extrapolation of the X-ray power law. 

The second part of the calculation involves using the 
CLOUDY model in conjuction with the observed blob line 

luminosity to calculate A c ioud, and from this the model- 

(2) 

predicted rate, N io ^, at which hydrogen-ionizing photons 

are absorbed. For the case of ionization-bounded clouds, the 

(2) 

^ion thus calculated is constant, regardless of the assumed 
values of gas density and fscatt in the CLOUDY model; this 
follows from the Zanstra method (Osterbrock 1989) which 
dictates that each absorbed hydrogen-ionizing photon leads 
to the production of one Balmer line photon, and the H/3 
luminosity of the blob is known observationally. Thus for 
ionization-bounded clouds, setting TV' 2 ' = leads to the 
relation A c i ou d = 36/ fscatt- The projected area of the blob 
as we define it is ~ 18 kpc 2 , so on the assumption that it 
has at most a similar extent along the line of sight and is 
probably composed of individual cloudlets with a covering 
factor of less than unity, we take A c i ou d < 18, which implies 
fscatt > 2. Given that Ogasaka et al. (1997) deduce that 
the small observed value of Lx/LboI suggests fscatt ~ 10, 
the conclusion to be drawn from this calculation is that we 
cannot rule out that the blob is ionized by continuum radia- 
tion from the active nucleus. [Models with matter-bounded 
clouds do not absorb all of the incident ionizing continuum, 
and lead to larger values of A c ioud for a given fscatt, and 
vice versa] . 



3.2 The ISIS data 

A high resolution long-slit spectrum of this object was ob- 
tained for a number of reasons. Firstly, and as mentioned 
at the end of section 3.1.1., the parameters extracted from 
the fits to the [OIII]+H/3 complex of the ARGUS data could 
be subject to a host of systematic effects. Secondly, the in- 
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Figure 6. The nuclear spectrum of IRAS F20460+1925 
taken with the ISIS spectrograph on the WHT. 



ability to decompose the [NII]+Ha profile prohibits the de- 
termination of any intrinsic extinction. Finally, the deterio- 
ration of the signal shortward of H/3 prevents the measure- 
ment of additional line ratios which could aid in the dis- 
crimination between shock+precursor ionization and AGN- 
photoionization. The ISIS long slit was oriented at a position 
angle of 30 degrees, in order to include a knot of emission 
which was noticed in an archival HST WFPC2 V-band im- 
age of the object (PID 5463). At a projected nuclear distance 
of 30 kpc this feature could be the tidal remnant of a merger. 

The nuclear ISIS spectrum is shown in Fig. |i| the two 
velocity components to the [OIIIJ+H/3 profile are clearly re- 
solved, with the blueshifted one being the more intense and 
the more highly ionized. A two component fit yields veloc- 
ity widths of 870 ± 40 and 1010 ± 20 kms" 1 FWHM for 
the systemic and blueshifted components respectively, and 
a separation of ~ 990 kms -1 between them. The similarity 
in these widths means that the parameters obtained from 
the common-width model fits to the ARGUS data will not 
suffer seriously from any of the systematic effects described 
in Appendix A. This is borne out by the fact the velocity off- 
sets of the blueshifted components derived from the ARGUS 
and ISIS data are in near-agreement, as are the [OIII]/H/3 
and [OIII]/[OII] ratios for the two components. The latter 
also gives independent confirmation of the ARGUS flux cal- 
ibration blueward of H/3. 

The spatial variation in the intensity of the blueshifted 
component along the direction of the slit is consistent with 
its being spatially unresolved, thus confirming what was sug- 
gested by the ARGUS data. Given that the new data were 
taken under seeing of ~ 1 arcsec, it can be deduced that the 
true surface brightness profile of the blueshifted component 
is peaked within a radius of 0.5 arcsec [= 2 kpc] 

The internal extinction was estimated using two- 
component fits to the H7 emission line. The velocity widths 
and separation of the two components were fixed at the val- 
ues obtained from the [OIIIJ+H/3 fits. When compared with 
the theoretical case B value for an electron temperature of 
10 4 K, the H/3/H7 value is consistent with an extinction of 
Ay ~ 1.8 mag for both the systemic and blueshifted com- 
ponents in the nuclear spectrum. Off nucleus, the systemic 



H/3/H7 is, within the errors, consistent with case B and 
therefore with no extinction, whilst that of the blue com- 
ponent remains close to the nuclear value (consistent with 
its being unresolved). For comparison, Frogel et al. (1989) 
deduced an extinction of Ay ~ 2.4 mag from the opti- 
cal narrow lines, and Veilleux et al. (1997) used their Paa 
flux (broad+narrow line) and the narrow line Ha flux to 
compute an upper limit to the extinction to the NLR of 
Ay — 4.6 mag. In the diagnostic diagram of Fig. || we show 
the correction factor which must be applied to the observed 
[OII]/[OIII] from the ARGUS data in order to correct for 
Ay — 2 mag. 

In addition to the Balmer lines and those of [Oil] 
and [OIII], the following emission lines could be dis- 
cerned in the nuclear spectrum: HeIIA4686, [SII]A4069,4076, 
[NeIII]A3965,3976 blended with He, [NeIII]A3869,3889 
blended with H£, and [NeV]A3426. Because of blending and 
the presence of multiple velocity components, velocity de- 
composition was attempted for Hell and [NeV] only. The 
position of the former is consistent with its being pro- 
duced predominantly by the blueshifted component and 
where it could be detected (within ~ 0.5 arcsec of the nu- 
cleus) the HeII/H/3 ratios were found to be 0.24 ± 0.08 and 
0.14 ±0.05 for the blueshifted and systemic components re- 
spectively. The blueshifted component thus lies at the same 
position (in terms of the values of the shock velocity and 
magnetic parameter) in both Fig. ^ and the HeII/H/3 ver- 
sus [OIII] /H/3 diagnostic diagram of Dopita & Sutherland 
(1995) for shock+precursor ionization. This constitutes self- 
consistent evidence for a shock+precursor ionization mech- 
anism for the blueshifted component, although the implied 
model shock velocity of around 500 kms" 1 is at variance 
with the observed kinematics. The [NeV]A3426 line could 
be resolved into systemic and blueshifted components, hav- 
ing [NeV] /H/3 values of 0.085 ± 0.03 and 0.28 ± 0.1, respec- 
tively. When an upward extinction correction factor of 2 
is applied, the value for the blueshifted component is just 
compatible with the value of 0.39 predicted by the highest 
velocity (500 kms" 1 ) shock+precursor model of Dopita & 
Sutherland (1995) for some fiducial value of the magnetic pa- 
rameter (they do not present a diagnostic diagram involving 
[NeV]). 

The [OIII]/H/3, [OII]/[OIII] and HeII/H/3 values for the 
systemic component on nucleus satisfactorily match the pre- 
viously discussed CLOUDY models for ionization parame- 
ters U ~ 10" 3 . The observed strength of [NeV]A3426 is, how- 
ever, much stronger than these models predict, but this is a 
well known problem with the nuclear spectra of active galax- 
ies (Binette, Courvoisier & Robinson 1987; Binette, Robin- 
son & Courvoisier 1987). We note that it can be solved by 
appealing to continua with additional high-energy ionizing 
radiation (above the Ne 3+ ionization potential of 97 eV), or 
to the presence of some matter-bounded clouds which pro- 
duce stronger lines from the more ionized species. 

Finally in connection with this object, we mention that 
no emission (neither line nor continuum) was seen at the 
position of the knot noticed in the HST image and with 
which the slit was aligned. 
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3.3 Summary of the ARGUS and ISIS data on 
IRAS F20460+1925 

The ARGUS data revealed that the gas kinematics in this 
object are violent and complex. These findings were con- 
firmed by the higher resolution ISIS data which showed that 
there are two velocity components to the [OIII]+H/3 complex 
in the nuclear region, with FWHM of 870 ±40 and 1010 ±20 
kms -1 , separated by ~ 990 kms' 1 . The blueshifted of the 
two appears to be spatially unresolved and could be due to 
an outflow of some kind, perhaps similar to the starburst- 
driven 'superwinds' seen in some lower luminosity IRAS 
galaxies or to the minor axis outflows seen in at least a 
quarter of Seyfert galaxies (Colbert et al. 1996). With the 
exception of two areas of redshifted emission to the west 
and north-east of the nucleus, the systemic velocity field is 
otherwise relatively devoid of structure. 

The ionization state of gas at the systemic veloc- 
ity is more consistent with photoionization by nuclear 
continuum radiation than by shocks and their precur- 
sors. The blueshifted component is more highly ionized 
and its emission-line intensity ratios are consistent with 
shock±precursor ionization, although the implied shock ve- 
locity cannot easily be reconciled with the observed kine- 
matics. The north-eastern blob also stands out as an area of 
high excitation in maps of [OIII]/[OII] and [OIII]/H/3 and 
simple modelling suggests that it could be ionized by con- 
tinuum radiation from the active nucleus; it may therefore 
resemble the photoionization cones seen in some lower lumi- 
nosity Seyfert galaxies. 



4 RESULTS FOR IRAS F23060+0505 

IRAS F23060±0505 is the second QSO-like object in the 
sample. At a redshift of 0.174, its inferred 0.5-120/^m spec- 
tral energy distribution and luminosity (2.9 x 10 46 ergs -1 ) 
are characteristic of a Seyfert 1 galaxy, but its optical emis- 
sion line spectrum more closely resembles that of a LINER 
or Seyfert 2 (Hill et al. 1987; Hough et al. 1991). Various 
authors have utilised near-IR/optical spectroscopy and po- 
larimetry to model its properties as a hidden Seyfert 1 nu- 
cleus: Hough et al. (1991) accounted for the wavelength de- 
pendence of the polarization with a dichroic model in which 
the nucleus and broad-line region (BLR) are viewed directly 
through aligned dust grains. Subsequent spectropolarimetry 
by Young et al. (1996b) led to a replacement model in which 
light from the BLR is scattered towards us by electrons in 
a cone with a half-opening angle of 45 degrees, inclined at 
50 degrees to the line of sight. The scattered radiation un- 
dergoes extinction equivalent to Av = 2.9 mag, with ex- 
cess flux in the near-IR accounted for by a dichroic view 
to the near-IR-emitting region through Av = 20 mag. The 
direct view to the BLR is obscured by Av = 30 ± 4 mag. 
Hines (1991) and Veilleux et al. (1997) detected broad Paa 
with FWHM ~ 2900 kms -1 , and derived lower limits on the 
visual extinction to the BLR of Av > 4.4 and > 7.4 mag, 
respectively. 

The findings of Young et al. (1996b) were supported 
by the ASCA 2-10 keV spectroscopy of IRAS F23060±0505 
reported by Brandt et al. (1997). The latter authors demon- 
strated that the data were fitted best by a 'two light- 
path model' in which the received X-rays comprise those 



coming directly to us through a column density of A^h — 
8 x 10 22 cm -2 and an additional electron-scattered compo- 
nent (itself absorbed by Nn — 5 x 10 21 cm -2 ). The covering 
fraction of the obscuring torus (as seen from the central 
source) was deduced to be 99.2 per cent. The unabsorbed 
spectrum is a power-law with a photon index of 2 and there 
is evidence for the presence of an iron K-line, both of which 
betray the presence of an AGN. 

The host galaxy was imaged by Hutchings & Neff (1988) 
and appears to be moderately blue, indicative of a relatively 
young stellar population. They suggest that the presence 
of a tidal tail composed mainly of old red stars attests to 
its being in the early stages of tidal disruption. Their VLA 
snapshot revealed that the system contains an unresolved 
radio source with a 20cm flux of 6.7 ± 0.3 mjy. In contrast, 
Heisler et al. (1994) noted that the source is unresolved at 
6cm but resolved at 20cm with fluxes of 3.0 and 13.5mJy, 
respectively (but they do not comment on the nature of the 
resolved structure). 



4.1 The ARGUS data 

4-1.1 The line- emitting gas: its kinematics and 
morphology 

Inspection of the CCD frames showed narrow line profiles 
asymmetrically skewed towards the blue, as reported by 
Young et al. (1996b) [cf. IRAS F20460±1925, with its two 
distinct, relatively symmetrical, components]. There is thus 
less physical justification for the fitting of a two component 
model for this object than there was for IRAS F20460±1925 
(still less, one in which both components have the same 
width). In the nucleus, the second component is blueshifted 
by ~ 830 kms" 1 with respect to the systemic and has ~ 20 
per cent of the latter's flux in [OIIIJA5007. Fig. ^ shows the 
[OIII]±H/3 complex in both the nucleus and a secondary 
enhancement in the blueshifted component (hereafter the 
north-western blob). 

Fig. |i] shows reconstructed images of the object in 
various continuum-subtracted emission lines. The surface 
brightness distribution of the systemic [OIIIJA5007 is cen- 
trally peaked with an approximately circular morphology, 
and there are two separate areas of blueshifted [OIII] emis- 
sion. The first peaks on the same (nuclear) fibre as the 
systemic with a surface-brightness profile which follows the 
instrumental PSF and is thus unresolved (this fitted com- 
ponent may not, however, be a physically distinct veloc- 
ity component). The second region of blueshifted compo- 
nent emission is offset from the nucleus by 7.5 kpc to the 
north-west. The total [OIII]A5007 luminosity of the system 
is 7.6 x 10 42 ergs -1 , 18 per cent of which is contributed by 
the blueshifted component. 

The [OII]A3727 fluxes were derived from fits to this 
line with single-component gaussians. The decomposition of 
the Ha±[NII] complex could be performed and was mod- 
elled with single component narrow lines and broad Ha, the 
width of the latter being fixed at the value 4550 kms -1 ob- 
tained by Young et al. (1996b), owing to the contamination 
of its blue-wing by atmospheric absorption. The morphol- 
ogy of the narrow-line Ha emission map differs from that of 
[OIII]A5007, the former being significantly more extended 
along a north-south axis. 
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Figure 7. The [OIHJ+H/3 complex in the nuclear fibre and 
the north-western blob in IRAS F23060+1925. The dashed 
line indicates the position of the systemic component in the 
nuclear fibre. 



See separate figure 



Figure 8. Reconstructed images of IRAS F23060+0505 
in the following emission lines: systemic and blueshifted 
components of [OIII]A5007, narrow and broad Ha, and 
[OII]A3727 (see text for details). The square root of the 
integrated line flux has been plotted for each fibre, in or- 
der to bring out more clearly the secondary enhancement in 
the blueshifted component of [OIII] 7.5 kpc to the north- 
west of the nucleus (the north-west blob, as indicated). The 
scale is in units of 10 -8 erg°' 5 cm -1 s -0 ' 5 . Discs are drawn 
only for those fibres with a significant line detection, and 
the actual flux ranges spanned by the plotted fibres are 
(in units of 10 -15 erg cm -2 s" 1 ): 0.0033-4.1 (for systemic 
[OIII]), 0.0026-0.87 (for blueshifted [OIII]), 0.03-0.77 (for 
[Oil]), 0.13-2.8 (for broad Ha) and 0.013-2.2 (for narrow 
Ha). North is at the top and east to the left. 



Figure 9. The upper map shows the FWHM common to the 
systemic and blueshifted components in a two-component 
fit to the [OIII]+H/3 complex of IRAS F23060+0505; the 
linewidth was unresolved in those fibres marked with as- 
terisks (see text). The lower map shows the radial velocity 
of the systemic component of [OIHJ+H/3 emission relative 
to the nuclear fibre (marked with a cross). The dashed line 
denotes the axis along which the velocity cut of Fig. [u] is 
taken. North is at the top and east to the left. 

Fig. shows the velocity of the systemic component 
of the [OIII]A5007 emission relative to the nucleus, and its 
fitted FWHM. Dipolar structure is apparent in the systemic 
velocity field, and that it resembles normal galactic rotation 
is evident from Fig. [ii^, in which we show a velocity cut along 
the axis indicated by the dashed line of Fig. ^. The rotation 
curve is asymmetric, in the sense that the kinematic centre of 
the galaxy is not coincident with the nucleus (defined as the 
peak in the [OIII]A5007 emission). This will be returned to 
in section 4.2. The region of increased line-width (FWHM^ 
600 km s -1 , as compared with ~ 400 km s -1 for the rest of 
the nebula) lies within, but is not exactly co-spatial with, 
the blueshifted half of the dipolar structure. This suggests 
that the increased line width is not a correlated response of 
our two-component model to the dipolar velocity structure 
(or vice- versa). 

In the interpretation of the fitted velocity field, we refer 
to the discussion of the profile-fitting systematics presented 
in Appendix A. The simulations described therein were per- 
formed in order to investigate how the parameters extracted 
from fits with a two component, common-width model de- 
pend upon the true functional form of the profile, in which 
the widths of the two components may be quite different 
from each other. The results confirm that the dipolar struc- 
ture in the systemic velocity is real (indeed, it could be seen 
on the raw CCD frame). In contrast, the fitted relative in- 
tensity and velocity offset poorly represent the properties of 
any second component. That the increased FWHM might in 
part be a fitting artefact (ie. not entirely due to a genuine 
increase in the FWHM of the systemic) was suggested by 
its coincidence with a contigious group of fibres in which 
the model significantly underpredicts the flux observed be- 
tween the [OIIIJA5007 and A4959 lines. The correlation of 
these two features was seen during the simulations to result 
separately from the following changes to the true functional 
form of the profiles (assuming that the blue wing to the line 
is produced by a second gaussian component, much broader 
than the systemic): 
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Figure 10. Radial velocity of the systemic component 
of [OIII] with respect to the nuclear fibre of IRAS 
F23060+0505, along the dashed line of Fig. § 
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Figure 11. Narrow line Balmer decrement plotted against 
the fitted line width for groups of fibres in IRAS 
F23060+0505. The arrow denotes a three sigma lower limit 
for a group of fibres in which H/3 is not detected 



(a) A fibre-to-fibre variation in the FWHM of the sec- 
ond component. This would imply that the latter is spa- 
tially resolved, but would not contradict the finding that the 
second fitted component is spatially unresolved because its 
fitted relative intensity significantly underpredicts the true 
value. It would, however, be hard to understand how a spa- 
tially resolved outflow could acquire a sufficiently large width 
(FWHM~ 1350 km s^ 1 ) for this effect to be observable. An 
unresolved outflow naturally generates a broad profile by 
integration over the spatial extent of the outflow with its 
spread in line of sight velocity. 

(b) A spatial variation in the relative intensities of the 
two components. This could arise if the blueshifted compo- 
nent were due to an unresolved nuclear outflow, superim- 
posed on the spatially-resolved systemic. Fibre-to- fibre vari- 
ations in the systemic flux could be intrinsic to the line- 
emitting gas or produced by preferential dust-obscuration. 
Indeed, an archival HST WFPC2 image of the object (PID 
5463) shows a prominent dust lane to the east of the nucleus, 
coincident with the region of interest. That this dust is in the 
form of a foreground screen or mixed with the line-emitting 
gas is confirmed by the narrow- line Ha/H/3 ratio. Fig. |ll| 
shows the Balmer decrement plotted against the fitted a for 
groups of fibres spread across the nebula: a clear correlation 
is apparent, suggesting that this systematic effect may be 
operating. 

Thus, it cannot be ruled out that the observed increase 
in the fitted FWHM is a systematic effect, as described in 
(b) above, and not entirely due to a genuine increase in the 
velocity dispersion of either the systemic or blueshifted com- 
ponent. This conclusion is insensitive to whether or not the 
second component is spatially resolved, which in any case 
is difficult to ascertain from the present data. We empha- 
sise, however, that the correlation in Fig. [y] may indicate 
the genuine co-existence of a higher velocity dispersion with 



increased extinction. We shall return to this issue in the 
discussion of the ISIS data in section 4.2. 

We discuss here the secondary enhancement in the 
blueshifted component to the [OIII] at the north-western 
blob. At this position it is the stronger of the two veloc- 
ity components, with a flux in [OIII]A5007 of (4.6 ± 1.3) x 
1CP 17 erg cm -2 s _1 . It is blueshifted relative to the nuclear 
systemic emission by 400 ± SOkms" 1 and a 3<r lower limit 
on [OIII] /H/3 is 4.7. In view of the poor signal-to- noise and 
doubts about whether the blueshifted component centred 
on the nucleus is spatially resolved or not, it is difficult to 
ascertain whether this feature is a discrete 'blob' or merely 
a region of enhanced surface-brightness in a more extensive 
blueshifted component. The kinematics are not helpful in 
this respect either, due to the uncertainty in the fitted ve- 
locity offset of the second component in fibres beyond the nu- 
cleus and the complications introduced by the fitting system- 
atics. In the HST V-band image the position of this 'blob' is 
seen to coincide with a local deficit in the surface-brightness. 
The morphology of the gas/stars around this 'hole' suggests 
that some kind of disturbance has recently occurred there. 

4-1-2 Ionization structure 

Maps depicting the spatial variation in the line ratios 
[OIII] /[Oil] and [OIII] /H/3 across IRAS F23060+0505 are 
presented in Fig. [H| The [OII]A3727 flux is that derived 
from a single gaussian fit to the line. The [OIII]A5007 flux 
for this map is taken to be the systemic component alone, 
because if, as is likely, the blue wing is more highly ionized 
than the systemic, its contribution to [Oil] will be much less 
than its ~ 20 per cent relative strength at [OIII]. 

The map of the systemic [OIII] /H/3 suggests some 
anisotropy in the excitation of the gas, with the ratio being 
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Figure 12. The upper figure shows [OIII] A5007/[OII] A3727 
for all fibres with detectable [Oil] in IRAS F23060+0505, see 
text for details. The lower figure shows [OIII]/H/3 (systemic 
component), with the outlined fibres being the only indi- 
vidual fibres in with detectable H/3; elsewhere fits were also 
carried out to groups of fibres (not outlined), with diamonds 
indicating 3cr lower limits. 



highest west and north-west of the nucleus, in the vicinity 
of the north-western blob. If this gas is photoionized by the 
nuclear continuum source, this could be due to the com- 
bined effects of an ionization cone (produced by collimation 
of the radiation field within the nuclear regions) and a gen- 
eral patchiness in the distribution of obscuring material in 
the galaxy itself (as the fact that the apex of the 'cone' is 
not centred on the nucleus would imply). The latter may 
also account for the offset between the kinematic centre of 
the galaxy and the observed position of the nucleus. 

In Fig. |l^ we show a line ratio diagnostic diagram for 
the individual fibres in Fig. [l^. The role of nuclear photoion- 
ization is again difficult to assess because of our ignorance 
of the shape of the optical- UV continuum. Of the two model 
loci shown in Fig. that corresponding to shock+precursor 
ionization provides the closer description of the ionization 
state. In the context of this model, fibres at greater radii 
are consistent with lower shock speeds. It is difficult to 
test shock+precursor ionization by looking for the expected 
power-law dependence of the H/3 flux on the shock-velocity 
because systematic effects in the spectral fitting and projec- 
tion effects undoubtedly mean that the fitted FWHM in a 
particular fibre does not directly indicate the speed of the 
shock responsible for the line emission. Nevertheless, shock 
excitation remains attractive, not least because X-ray and 
spectropolarimetric studies indicate that the nuclear contin- 
uum source is almost completely covered with thick obscur- 
ing material. As in section 3.1, we mention, however, that a 
single diagnostic diagram cannot be used to rigorously dis- 
criminate between the two ionization mechanisms. 



4.2 The ISIS data 

In acquiring an ISIS spectrum of this object we sought a 
better understanding of the origin of the increased FWHM 
and of the nature of the blue-asymmetry in the line profiles, 
as well as a greater number of emission-line diagnostics. A 
slit position angle of 135 degrees was chosen, in order to 
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Figure 13. [OIII]/H/3 versus [OIII]/ [Oil] for the systemic 
emission components of the fibres in IRAS F23060+0505 
where all the relevant lines are detected. The photoioniza- 
tion and shock+precursor loci and the extinction-correction 
arrow are as in Fig. [B|. Star symbols denote the seven central 
nuclear fibres, open circles fibres at a radius of 0.8 arcsec, 
open squares fibres at 1.2 arcsec radius and filled circles some 
fibres at radii > 1.5 arcsec. 
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Figure 14. The nuclear ISIS spectrum of IRAS 
F23060+0505. 



cover both halves of the dipolar velocity field and the north- 
western blob. 

Fig.jli] shows the nuclear ISIS spectrum, with its blue- 
asymmetric line profiles. Spectra along the slit were fitted 
with the two component model of the [OIII]+H/3 complex 
(with the two widths no longer constrained to be equal) and 
the variations in the fitted parameters are shown in Figs. |l5[ 
and 0. The principal findings are that: 

(1) The asymmetry in the systemic velocity rotation 
curve seen in Fig. ^ is confirmed. 
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Figure 15. Variation of the fitted model parameters with 
position along the ISIS slit for the systemic component in 
IRAS F23060+I925. From top to bottom are shown the 
[OIII]A5007 flux, the velocity with respect to the nucleus 
and the FWHM. The north-western and south-eastern ends 
of the slit are indicated. 



(2) The width of the systemic is constant on the SE 
half of the dipolar structure, and declines to a much smaller 
value on the NW half. 

(3) The velocity offset of the blueshifted component 
from the nuclear systemic is constant to the SE of the nu- 
cleus and decreases in amplitude towards to the NW. 

(4) The width of the blueshifted component decreases 
from SE to NW. 

(5) The [OII]/[OIII] ratio is constant to the SE of the 
nucleus, but increases NW of it. 

We thus conclude that, if the blue-asymmetric line pro- 
files are due to a single second component - and for the pur- 
pose of ascertaining whether this is the case the ISIS data are 
no more useful than the ARGUS data - it is poorly modelled 
by the common- width two-component model used to fit the 
latter. The revision of some of the conclusions of section 4.1 
is thus necessary. Firstly, it can no longer be maintained that 
any second component is spatially unresolved. Secondly, the 
systematic effects which were cited as possible causes for 
the region of increased FWHM in Fig. ^| are probably not 
primarily responsible for it, the ISIS data showing instead 
that there are genuine increases in the widths of both the 
systemic and blueshifted components. 

We discuss now the ionization state of the g 
revealed by the ISIS data. The measured [OIII]/H/3 and 
[OII]/[OIII] values for the systemic component accord well 
with the ARGUS values, despite the systematic effects in 
the fits to the latter. Some information on the ionization 
state of the blueshifted component, which could not be re- 
liably measured using the ARGUS data, is presented in 
Fig. |l7]; as the presence of a second component to H/3 is 
questionable, we show 3 sigma lower limits in the [OIII]/H/3 
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Figure 16. As in figure |l5| but for the blueshifted component. 
The velocity is relative to the nuclear systemic. 
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Figure 17. From top to bottom are shown the [OIII]/H/3 
ratios of the systemic and blueshifted components and 
[OII]/[OIII] for the blueshifted component, with arrows de- 
noting three sigma lower limits for spectra with a question- 
able second component to H/3. 



plot. The [OII]/[OIII] ratio is approximately constant in the 
blueshifted half of the dipolar velocity field, but increases 
with nuclear distance on the other half of this structure. 

Whilst many of the other lines present in the IRAS 
F20460+1925 spectrum also appear in that of IRAS 
F23060+0505, the analysis of them is here less straight- 
forward. This is because the signal-to-noise ratio is lower, 
and the disparity in their widths means that the blueshifted 
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wing falls below the noise level before the systemic part of 
the line. The smaller separation of the two components also 
makes it difficult to attribute unambiguously a line which 
can be fitted with only a single component to one or the 
other of the two. A two component fit to HeIIA4686 shows 
that the blueshifted component is marginally stronger, and 
yields HeII/H/3 ~ 0.1 - 0.23 and (lower limits of) ~ 0.3 - 
0.9, respectively, for the systemic and blueshifted compo- 
nents. Similarly, a fit to [NeV]A3426 is dominated by the 
blueshifted component and yields [NeV]/H/3 — 0.1 — 0.18 
and (lower limits of) ~ 0.5 — 1.1, respectively, for the 
two components. The systemic HeII/H/3 and [NeV]/H/3 val- 
ues are, when the latter is extinction-corrected, consistent 
with the shock+precursor models of Dopita & Sutherland 
(1995) (as the ARGUS data suggested), with velocities up 
to ~ 500 kms" 1 . The latter range matches that shown in 
Fig. [t| for the systemic FWHM. 

The observed line-ratios of the blueshifted component 
are less-easily understood; the measured [OII]/[OIII] and 
[OIII]/R3 ratios do not discriminate between the two ion- 
ization mechanisms and the HeII/[OIII] and [NeV]/[OII] ra- 
tios fail to match the predictions of the otherwise plausible 
CLOUDY models or the values given by Dopita & Suther- 
land (1995) for the highest velocity shock+precursor models. 
Given, however, that the latter is for 500 kms" 1 whereas 
the FWHM of the blueshifted component is much larger 
(Fig. ^), it seems more likely that the gas is ionized by 
high velocity shocks (with velocities ~ 1000 kms" ) than 
by the non-thermal continuum of a central active nucleus. 

The H/3/H7 ratio for the blueshifted component is 
poorly constrained, and that of the systemic is consistent 
with an extinction of Ay — 1.2 mag on nucleus. Fits to 
average spectra in the red- and blueshifted halves of the 
dipolar structure yield extinctions of Ay ~ 0.4 and 1.1, re- 
spectively, although the extinction on the blueshifted half 
is not uniform, the Balmer decrement being in places con- 
sistent with the case B value. Nevertheless, that the ex- 
tinction is generally higher on the blueshifted half of the 
structure is consistent with the findings of the ARGUS data 
and with the location of the dust lane in the HST image 
of this object. For comparison, Hough et al. (1991) deduced 
Ay = 3.4±0.5 mag from single component fits to the narrow 
Ha and H/3 lines, and Veilleux et al. (1997) used their com- 
bined broad+narrow Paa flux to estimate an upper limit on 
the NLR extinction of Ay < 5.5 mag. 

Although the slit was aligned to pass through the sec- 
ondary enhancement in the blueshifted component noticed 
in the ARGUS data, the noise level in the ISIS data is such 
that nothing more can be added to the discussion of sec- 
tion 4.1.1. 

4.3 A kinematic model for IRAS F23060+0505 

In the context of the merger-induced formation scenario, we 
speculate here on an interpretation of this object's kinemat- 
ics. The long tidal tail (Hutchings & Neff 1988) signals an 
on-going merger in which the galaxies are near to the final 
coalescence. If not the result of dust obscuration, the asym- 
metric rotation curve may be due to the merger-induced 
disturbance of the pre-encounter velocity field of the princi- 
pal progenitor galaxy; the increased FWHM and dust lane 
in the same region may have a similar origin. The kine- 



matics of the blueshifted component shown in Fig. |l6| are 
suggestive of some sort of bipolar nuclear outflow. Within 
the ULIRG evolutionary paradigm, the latter would be as- 
cribed to a starburst-driven superwind, but we also note that 
the spatial structure of the [OIII] profiles may be related to 
the interaction of the previously-noted radio structure with 
the ambient medium, similar to that found by Whittle et 
al. (1988) for a number of Seyfert galaxies (see section 7.2). 

4.4 Summary of the ARGUS and ISIS data on 
IRAS F23060+0505 

The systemic velocity field of this object has a dipolar struc- 
ture with an axis along a position angle of 30 degrees. The 
rotation curve is, however, asymmetric, in the sense that 
it has amplitudes of 100 and 200 kms -1 on its blueshifted 
and redshifted halves, respectively. The systemic FWHM is 
approximately constant at 475 kms" 1 in the SE half, but 
declines towards the NW. The second velocity component 
with which the blue asymmetry in the lines of the [OIII]+H/3 
complex was modelled is much broader than the systemic 
(FWHM up to 1200 kms" 1 ) and offset from the latter by 
-600 kms" 1 on nucleus. A secondary enhancement to the 
blueshifted velocity component is seen 7.5 kpc NW of the 
nucleus, but the data do not allow us to ascertain whether 
it is a discrete 'blob' of emission or a region of enhanced 
surface-brightness in a more extensive blueshifted compo- 
nent. It would be of interest to obtain a high-resolution ra- 
dio map of this object, to examine the relationship between 
the spatial structure of the [OIII] profiles and any extended 
radio emission. 

The [OIII] /H/3 and [OIII] /[Oil] values from the AR- 
GUS data suggest that shocks and their precursors ionize 
the systemic velocity component. Other diagnostics from 
the ISIS data confirm this, and show that the implied shock 
velocities match the observed systemic FWHM. The ion- 
ization state of the blueshifted component is not consis- 
tent with the CLOUDY model used throughout this paper, 
nor with the highest velocity (500 kms -1 ) shock+precursor 
model of Dopita & Sutherland (1995). Since the blueshifted 
FWHM is much greater than 500 kms" 1 and studies at 
other wavelengths indicate that the active nucleus is very 
heavily obscured, we also favour shocks (perhaps faster than 
those which Dopita & Sutherland compute) as the ionization 
mechanism for the blueshifted component. 



5 JKT MULTI-BAND IMAGES OF IRAS 
F20460+1925 AND F23060+0505 

The ARGUS and ISIS spectroscopic data have enabled us 
to study in some detail the kinematics and ionization struc- 
ture of the ionized gas in these systems. Further clues to 
their evolutionary history are provided by the colours of the 
host galaxy stellar populations. Hutchings and Neff (1988) 
compared the B- and R-band images of IRAS F23060+0505 
and observed a series of arms or curved luminous regions 
extending in all directions from the nucleus, which is itself 
resolved at radii of < 1 arcsec. They observed also a series 
of knots or secondary nuclei, particularly in B light, at dis- 
tances of 2 to 4 arcsec from the nucleus, and claimed that 
these colour differences distinguish regions of recent massive 



© 0000 RAS, MNRAS 000, 000-000 



16 



star-formation from older stars. For IRAS F20460+1925, 
Frogel et al. (1989) obtained images in the V- and Fbands 
but did not report any colour gradients across the galaxy, 
which has an absorption-corrected Mv = —22.3, compara- 
ble to that for giant spiral or elliptical galaxies. They did, 
however, draw attention to the knot of V-band emission 6.5 
arcsec to the south-west of the nucleus (seen also in the 
HST V-band image, as mentioned in section 3.1) and de- 
duced that it, or any other companion object, must be at 
least 3 magnitudes fainter in I than the IRAS source. 

Using the JKT data, we confirm the complex colour 
structure of IRAS F23060+0505 reported by Hutchings & 
Neff, the most prominent features being the knot of B- 
band emission in the plume to the SW and a "fan-like" 
R-band extension of the nucleus in roughly the same direc- 
tion. The host galaxy of IRAS F20460+1925 is most appar- 
ent in the I-band and exhibits no discernable colour gradi- 
ents which, along with its quasi-elliptical surface-brightness 
profile, point to a relatively ancient stellar population. The 
V-band emission knot is absent from all three of the present 
filter images and from the ISIS spectrum, suggesting that its 
SED is peaked strongly in the V-band. Whilst these facts 
may be used as a basis for speculation about its proper- 
ties, it should be borne in mind that, lacking a knowledge 
of its redshift, its physical association with the main galaxy 
remains no more than a hypothesis. 



6 THE LOWER LUMINOSITY OBJECTS 

With bolometric luminosities close to 

10 13 L Q , IRAS F20460+1925 and F23060+0505 straddle 
the conventionally-adopted ultra/hyperluminous divide. In 
this section, we discuss the ARGUS observations of two less 
powerful objects [Lfir = 10 12 ' 2 and 10 12 ' 4 L©], included in 
the programme for comparison with the QSO-like systems. 
As mentioned in section 2.1, the seeing progressively deteri- 
orated during this period, thus precluding a comparison on 
the terms which had been originally envisaged. Whilst the 
study of the spatial structure of the kinematics and ioniza- 
tion state is slightly impaired, integrated properties such as 
the total line luminosities can still be extracted. 



6.1 IRAS F01217+0122 

IRAS F01217+0122 lies at z=0.137 and was included in the 
long-slit spectroscopic sample of Armus, Heckman & Miley 
(1989) (hereafter AHM), who commented upon its broad, 
blue-asymmetric emission line profiles. An R-band image 
shows two apparently interacting galaxies, but spectroscopy 
reveals that the larger of these is a foreground object (Ar- 
mus, Heckman & Miley 1990). 

Because of the poor seeing, the ARGUS spectra were 
analysed in hexagonal groups of seven fibres out to a ra- 
dius of 1.6 arcsec from the nucleus, beyond which no emis- 
sion line flux could be detected. The [OIII]+H/3 complex in 
these spectra were fitted with the two-component common- 
width model, with a series of four defective CCD columns 
on the blue- wing of the [OIII]A4959 line being masked out 
(but this does not adversely affect the fitted parameter val- 
ues). The kinematics appear to be almost entirely devoid of 



structure: the systemic velocity and fitted FWHM are sub- 
stantially constant across the nebula, the latter assuming a 
value of ~ 1050 kms" 1 . The velocity offset of the spatially- 
unresolved blueshifted component is ~ 1320 kms -1 . It may 
be the case that the observed kinematics are dominated 
by the effects of the seeing on a centrally-peaked surface- 
brightness profile, and that they are thus a poor represen- 
tation of the true velocity field. The total luminosity of the 
system in [OIII]A5007 is 3.53 x 10 42 ergs" 1 , 13 per cent of 
which comes from the blueshifted component. 

The nuclear spectrum (shown in Fig. ^) exhibits a 
number of strong emission lines blueward of H/3, includ- 
ing [OII]A3727, [NeIII]A3869, [OIII]A4363 and HeIIA4686 
which, when fitted with single component gaussian pro- 
files, yield line ratios of [OII]/[OIII]A5007 = 0.25 ± 0.03, 
[NeIII]/[OII]= 1.0 ± 0.16, HeII/H/3 = 0.53 ± 0.12 and 
[OIII]A4363/[OIII]A5007 = 0.019 ± 0.01 (the [OIII]A5007 
and H/9 fluxes being those of the systemic component). Note 
the unusually strong [NeIII]A3869 and [OIII]A4363; the lat- 
ter equates to an electron temperature of ~ 22000 K af- 
ter correction for an extinction of Av ~ 3.6 mag (de- 
duced from the Balmer decrement Fua/Fup = 9.46 given 
by AHM). The fitted value of [OIII]/H^ = 14.2 ± 2.3 on 
nucleus is consistent with the value measured by AHM, 
which was the largest of all the objects in their sample. 
Because there is evidence for underlying Balmer absorp- 
tion in some of the off-nuclear spectra of this object (in eg. 
Fig po| ), the nuclear spectrum may also suffer some contam- 
ination. When comparing the observed emission line inten- 
sity ratios with those predicted by models, we thus restrict 
our attention to those not involving H/3. The large value of 
[NeIII]A3869/[OII]A3727 strongly favours AGN photoioniza- 
tion over any of the shock+precursor of Dopita & Sutherland 
(1995); in particular, the CLOUDY model shown in Fig. ^ 
with a hydrogen nucleon density of 10°' 5 cm -3 and ioniza- 
tion parameter U ~ 10 -1 ' 5 , is able to reproduce the observed 
value of this ratio, whilst also being consistent with the mea- 
sured [OIII]A4363/HeIIA4686 and [OIII]A4363/[OIII]A5007 
values. The ionization state of the blueshifted component 
cannot be constrained at the signal-to-noise level of these 
data. 

That there may have been a relatively recent (< 1 Gyr), 
localised episode of star-formation is also suggested by the 
existence of a colour gradient across the galaxy. To quantify 
the strength of the blue continuum in excess of that from a 
giant elliptical galaxy (which, we note, may not be appro- 
priate for the underlying host galaxy) we fitted the spectra 
with a series of templates representing 05, B5, A5, F5 and 
G5 stellar models (Kurucz 1979). The spectra of the fibre 
groups were placed on a rest wavelength scale and changed 
to a dispersion of 15A per pixel to improve the signal-to- 
noise ratio, but were not corrected for intrinsic reddening, 
lacking as we do a reliable measure of the latter (because 
of the Balmer line absorption). Fitting was performed be- 
tween 3600-5400A (beyond which the template spectra are 
not tabulated) with strong emission lines being assigned zero 
statistical weight. As the smoothed galaxy spectra exhibit 
prominent H7 and H<5 absorption lines, the fitted models are 
not surprisingl y d ominated by stars of spectral types A and 
F. Figs, [it] and |20| compare the fitted spectra of two groups of 
off-nucleus fibres. The excess blue light of the entire galaxy 
(compared with an old elliptical galaxy stellar population) 
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Figure 18. A portion of the nuclear spectrum of IRAS 
F01217+0122. Note the strength of [NeIII]A3869. A series 
of bad columns on the blue- wing of the [OIII]A4959 line was 
masked out, as described in the text. 
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Figure 20. As in Fig. |l^ but for a group of fibres 1.2 arcsec 
north-west of the nucleus. The stellar model is dominated 
by A5 stars with a lesser contribution from those of class 
F5. 
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Figure 19. The spectrum of a group of fibres 1.2 arcsec south 
of the nucleus of IRAS F01217+0122, plotted at a dispersion 
of 15A per pixel on a rest- wavelength scale. The best-fit stel- 
lar template model is plotted up to 5400 A, and dominated 
by F5 stars. The strong emission lines were assigned zero 
statistical weight in the fit. 

can essentially be reproduced with ~ (1.8 ± 0.8) x 10 5 A5 
stars and ~ (1.4±1.1) x 10 10 F5 stars. The large errors, how- 
ever, prevent this result from being used to infer significant 
information about the star-formation history of the galaxy. 
Given that the main sequence lifetimes of 05, B5, A5 and F5 
stars are ~ 3.7 x 10 6 , 3.9 x 10 7 , 6.0 x 10 s and 3.3 x 10 9 yr, re- 
spectively (quoted by Shapiro & Teukolsky 1983) , it suffices 
to say that the results are consistent with an intermediate 
age (~ lGyr) stellar population. 

6.2 IRAS F01003-2238 

IRAS F01003-2238 lies at a redshift of 0.118 and has at- 
tracted the attention of numerous authors since its inclu- 
sion in the original sample of warm ULIRGs by Sanders et 
al. (1988b). Its popularity stems from its classification as 



a Wolf-Rayet galaxy by Armus, Heckman & Miley (1988) 
(hereafter AHMII). They interpreted a band of features at 
a rest wavelength of 4660A and other broad but weaker fea- 
tures as arising from ~ 10 5 late Wolf-Rayet stars of the W- 
N sub- type. More recently, Surace et al. (1998) presented 
WFPC2 B- and I-band images of this and nine other warm 
ULIRGs. The only notable circumnuclear feature they saw 
was a chain of knots extending from the SE to the NW 
within 1 arcsec of a putative nucleus. The latter are espe- 
cially prominent in the B-band and, if interpreted as regions 
of recent star-formation, have a mean age of about 7 Myr, 
consistent with the young ages ascribed to the stellar pop- 
ulations by AHMII. Despite being the brightest galaxy in 
a small group, it shows no discernable distorted structure. 
Near- infrared spectroscopy by Veilleux et al. (1997) found 
the galaxy to have a strong component of narrow Pact, with a 
blue-asymmetric line-profile base, but they saw no evidence 
for a hidden BLR. 

The gas kinematics were deduced from fits to the 
[OIII]+H/3 complex and are more globally structured than 
those of IRAS F01217+0122, showing a dipolar systemic 
velocity field with an amplitude of ~ 140 km s _1 across 
the galaxy. The FWHM of the lines are also broader 
in the blueshifted parts of this structure (as for IRAS 
F23060+0505), with a FWHM varying from ~ 450 to 
1000 km stover the same region. Both these features are 
shown in the kinematic maps of Fig ^l|. The second com- 
ponent to the [OIIIJ+H/3 complex used to model the blue 
wing to the lines is offset from the systemic by ~ —1150 
kms~ 1 on nucleus, with a relative intensity in [OIII] of 
~ 0.5. The luminosity of the entire system in [OIII]A5007 
is 1.32 x 10 42 ergs -1 , of which 29 per cent is contributed 
by the second component. 

The nuclear [OIIIJ/H/3 ratio of this object is the low- 
est in the present sample, with a value of 4.1 ± 0.8 The 
only other emission lines which can be discerned blueward 
of H/3 are [OII]A3727, and, with slightly less than 3cr signifi- 
cance, [Nelll] A3869. When fitted with single gaussians these 
yield line ratios of [NeIII]A3869/[OII]A3727 = 0.26 ± 0.15 
and [OIII]A5007/[OII]A3727 = 1.65 ± 0.6 (using the flux of 
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See separate figure 



Figure 21. The upper map shows the FWHM common to 
the systemic and blueshifted components in a two compo- 
nent fit to the [OIIIJ+H/3 complex of IRAS F01003-2238, 
corrected for the position-dependent instrumental resolu- 
tion. The fitting was carried out on hexagonal groups of 
seven fibres because of the poor seeing; the crosses denote 
the nuclear group and the asterisks mark a group in which 
the linewidth was unresolved. The lower map shows the ra- 
dial velocity of the systemic component of [OIIIJ+H/3 emis- 
sion. North is at the top and east to the left. 




4000 5000 6000 7000 8000 



Rest wavelength (A) 

Figure 22. The spectrum of the group of fibres in IRAS 
F01003-2238 approximately coincident with the star-clusters 
identified by Surace et al. (1998), plotted at a dispersion 
of 15A per pixel on a rest-wavelength scale. The best-fit 
stellar template spectrum, which is dominated by 05 stars, is 
plotted up to 5400A. The strong emission lines were assigned 
zero statistical weight in the fit. 

only the systemic [OIII] component). The fitting of a two 
component model to the [NII]+Ha complex yields a nuclear 
[Nil] /Ha ratio of 0.19±0.04, placing the object in the HII re- 
gion part of the [Nil] /Ha versus [OIII] A5007/H/3 diagnostic 
diagram of Veilleux & Osterbrock (1987). The low signal-to- 
noise of the present data precludes the identification of the 
Wolf-Rayet features noted by AHMII and the determination 
of the ionization state of the blueshifted velocity component. 

Due to the poor seeing (FWHMc^ 2 arcsec), the four 
knot-like star clusters and the putative nucleus identified 
by Surace et al. (1998) all fall within the seeing disk. It is 
thus not possible to use the ARGUS data to address the 



questions raised by the latter authors, concerning whether 
the observed Seyfert activity is associated with the puta- 
tive nucleus and whether the blue wing to the emission lines 
originates in a superwind from the star clusters. There is, 
however, a marked colour gradient across the galaxy, with 
the spectra of groups to the north and west of the nucleus (ie. 
in the vicinity of the star-clusters) possessing very steeply 
rising blue continua, clearly indicative of a young stellar pop- 
ulation. Stellar template spectra were used to quantify the 
amount of blue light using the procedure described above 
for IRAS F01217+0122, the only difference being that an 
extinction correction was applied to the IRAS F01003-2238 
spectra using the measured Ha/H/3 ratio. On nucleus, the 
observed ratio translates into Av = 1.4 mag, whilst on the 
star clusters it is consistent with the case B value and there- 
fore with no extinction. To the south and east of the nucleus 
extinctions as high as Av = 2.5 mag are deduced. This spa- 
tial gradient in the extinction conflicts with the findings of 
Surace et al. (1998), who claimed that the observed colour 
structure was indicative of a uniform dust screen. For com- 
parison, Veilleux et al. (1997) quote extinctions to the NLR 
derived from optical and near-infrared lines of Ay = 2.7 and 
5.5 mag, respectively. 

The spectra of the star clusters are dominated by 05 
and B5 stars (in terms of the fraction of the 4500A contin- 
uum flux they contribute), although all but the 05 contri- 
bution is poorly constrained in the fits. If fitted solely with 
05 stars, a total of 3.2 ± 0.2 x 10° such stars are required 
to account for these spectra. Fig. ^ shows the fit of the 
stellar template model to the spectrum of the two groups 
of fibres to the north and west of the nucleus, in the vicin- 
ity of the star-clusters. Considering that the main sequence 
lifetime of an 05 star is ~ 3.6 x 10 6 yr, this result is consis- 
tent with the age of 10 Myr deduced by Surace et al. (1998) 
for the oldest of these clusters by analysing evolutionary 
tracks on an M B / versus B-I diagram. When extinction cor- 
rections are applied to the spectra of the other groups of 
fibres, they too are seen to have similarly strong blue con- 
tinua. This suggests that star formation has occurred re- 
cently throughout the galaxy, and not just in the visible 
clusters where the extinction is observed to be consistent 
with Av = 0. On the assumption that the putative nucleus 
identified by Surace et al. (1998) is in fact a compact star- 
burst, we find that the extinction-corrected blue continuum 
of the entire galaxy can be accounted for with 9.2 ±0.5 x 10 6 
05 stars. Taking 10 5 ' 7 L© for the bolometric luminosity of an 
05 star (Allen 1973), this quantity of stars could plausibly 
account for the entire bolometric luminosity of the galaxy. 
This, along with the relatively low-ionization emission-line 
spectrum described above, may be taken to argue against 
the presence of an active nucleus in this object. Surace et 
al. (1998) found that the position of the putative nucleus on 
the M B ' versus B-I diagram was anomalous with respect to 
the other objects in their sample, and inconsistent with its 
being a reddened QSO nucleus. 



7 DISCUSSION 

Having dealt separately with each of the four ULIRGs in 
turn, we discuss here the comparative properties of the sam- 
ple as a whole. The findings of the preceeding sections are 
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Table 3. Summary of the properties of the four ULIRGS 
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★ See Table 2 

f BS refers to the velocity component blueshifted with respect to the systemic 
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summarised in Table 3. We aim to understand whether 
the observed diversity of kinematic and ionization struc- 
tures within our sample (not necessarily for the class of 
ULIRGs as a whole) can be accommodated within some gen- 
eral framework, such as the evolutionary scenario of Sanders 
ct al. (1988a). 



7.1 The systemic velocity fields 

Within our small sample of four objects, there appear 
to be two essentially different kinds of systemic velocity 
field. The first group comprises IRAS F20460+1925 and 
F01217+0122, in which the systemic velocity is approxi- 
mately constant across the galaxy with a uniformly high 
FWHM (~ 1000 kms -1 ); in contrast, objects in the sec- 
ond group (IRAS F23060+0505 and F01003-2238) exhibit 
more organised velocity fields, with dipolar structure and 
a FWHM which is smaller than in objects in the previous 



group but which is largest on the blueshifted half of the 
dipolar structure. In the context of the merger scenario for 
ULIRGs, we speculate that objects in the first group are old 
mergers in which the pre-encounter velocity structure has 
been destroyed by merger- induced turbulence, whilst those 
in the second are much younger and still retain some record 
of the velocity structure (eg. dominantly rotational) of the 
progenitors. It could instead be that the kinematic distur- 
bance is due a galactic-scale outflow (perhaps initiated by a 
merger) rather than to the effects of a merger per se. In a 
sample of this size, significance should probably not be at- 
tached to the fact that in both of the objects in the second 
group, the region of increased FWHM lies on the blueshifted 
half of the dipolar structure. 

In support of the above interpretation, we refer to the 
archival HST WFPC2 V-band images of IRAS F20460+1925 
and F23060+0505. That of the former shows approximately 
elliptical isophotes and, with the exception of a central point 
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source, an otherwise featureless surface brightness profile; a 
discrete knot of emission at a projected radius of 30 kpc may 
be a fossil tidal remnant. In contrast, the image of IRAS 
F23060+0505 shows a long tidal tail extending to the south- 
west, a prominent dust lane and numerous bright knots of 
emission, all of which suggest that the system is a rela- 
tively young merger which has recently undergone prodi- 
gious star formation, as claimed by Hutchings & Neff (1988). 
The JKT images (section 5) confirm these differences in the 
general morphology and colour structure of the two systems. 
Thus, IRAS F23060+0505 appears to be dynamically the 
younger of the two mergers (if indeed one did occur in IRAS 
F20460+1925), nominally consistent with the above expec- 
tations. 

For the lower-luminosity objects IRAS F01217+0122 
and F01003-2238 the above kinematic interpretation receives 
less direct support, as neither object shows obvious signs of 
interaction (Armus, Heckman & Miley 1990 and Sanders et 
al. 1988b, respectively). This is perhaps surprising, given 
that a very large proportion of ULIRGs reside in interacting 
or merging systems (Clements et al. 1996). Nevertheless, the 
large FWHM do suggest that mergers or galactic-scale kine- 
matic disturbances of some other kind (eg. superwinds) are 
occurring. As both the latter are intimately connected with 
bursts of star formation, we can compare the times elapsed 
since their onset by comparing the properties of their stel- 
lar populations. In section 6.1, the rising blue continuum 
of IRAS F01217+0122 was modelled with a ~ 1 Gyr old 
starburst, whereas that of F01003-2238 (section 6.2) was 
dominated by 05 and B5 stars, consistent with the age of 
10 Myr deduced by Surace et al. (1998) for the oldest star 
clusters. That IRAS F01003-2238 should have the younger 
stellar population thus conforms with the above interpreta- 
tion of the kinematics. 

7.2 The blueshifted velocity components 

The emission lines of the [OIII]+H/3 complex of all the 
objects are blue-asymmetric in appearance. We assumed 
that this is due to a physically distinct velocity compo- 
nent blueshifted with respect to the systemic part of the 
line (cf. the situation in which a single mechanism produces 
the entire line profile). In all but one of the objects (IRAS 
F23060+0505) it appears to be spatially unresolved, and 
(where reliable measurement is possible) more highly ion- 
ized than the systemic velocity component. Emission line 
intensity ratios point to ionization by shocks and their pre- 
cursors rather than to photoionization by an AGN con- 
tinuum. Table 3 shows that the [OIII]A5007 luminosity of 
the blueshifted component increases with LboI.fir (approxi- 
mately as LbJi.fir) an d that the outflow velocity on nucleus 
is generally higher in the less powerful systems. We caution, 
however, that these conclusions are potentially sensitive to 
the line-profile fitting systematics (Appendix A). 

Considering the paucity of constraints which we have 
upon their properties, we can do little more than speculate 
about the origin of these blueshifted components. One possi- 
bility is that they are outflows driven by the collective effects 
of many supernovae, the so-called starburst-driven super- 
winds; the latter are a common feature of lower- luminosity 
IRAS galaxies (see Armus, Heckman & Miley 1989, 1990) 
and are probably also present in other ULIRGs (eg. Veilleux 



et al. 1995; Kim et al. 1998). However, given the strong ev- 
idence for active nuclei in three of these objects, we should 
also consider their role. IRAS F20460+1925 is a case in 
point: unlike the other galaxies, the CCD frames of this 
object show a second, separate velocity component rather 
than a blue asymmetry to the line profile, which contributes 
a much greater fraction of the total [OIIIJA5007 than the 
blue wings of the other galaxies. And in addition, this sys- 
tem is likely to be an old merger with no sign of recent star 
formation to power a putative superwind. 

As mentioned earlier, Colbert et al. (1996) found that 
more than a quarter of Seyfert galaxies show evidence for 
large-scale minor axis galactic outflow and discussed the 
mechanisms which may drive them. Whittle et al. (1988) 
demonstrated that [OIII] profile substructure in Seyfert 
galaxies is frequently associated with extended radio struc- 
ture; such emission components typically have high excita- 
tion ([OIII]/H/3 ~ 10 — 15) and are undergoing systematic 
outflow from the galaxy. 

7.3 Ionization mechanisms 

Emission line intensity ratios suggest that in most of the ob- 
jects the gas at the systemic velocity is ionized by a hard ra- 
diation field (IRAS F01003-2238 is a possible exception - see 
below). Discriminating further between photoionization by 
an AGN and shock+precursor ionization is generally prob- 
lematic, as we are ignorant of the shape of the AGN contin- 
uum spectra and of the properties of the gas clouds. In IRAS 
F23060+0505 there is good evidence for shock+precursor 
ionization, from both the line ratios themselves and their 
concordance with the systemic FWHM derived from the ISIS 
data. As mentioned in section 4.1.2, this is consistent with 
indications from X-ray spectroscopy that the central engine 
is obscured over almost 47r steradian by a column density of 
./Vh ~ 10 23 cm~ 2 . The line of sight obscuration to the X-ray- 
emitting regions of IRAS F20460+1925 is also high, but less 
is known about its solid angle coverage; it is thus possibilc 
that radiation could escape to power the line emission, which 
has observed line ratios consistent with AGN photoioniza- 
tion. Furthermore, line ratio maps of this object show that 
the [OIII] /[Oil] ratio increases north-east of the nucleus, ter- 
minating in a region of high [OIII]/H/3. On account of its low 
FWHM, we speculated that this 'NE blob' could comprise 
gas in the host galaxy which is ionized by continuum radia- 
tion from the active nucleus, an interpretation which simple 
energetic considerations do not render implausible. These 
differences between IRAS F20460+1925 and F23060+0505 
may be taken to reinforce their previously-deduced relative 
positions along the merger-induced 'evolutionary sequence': 
the 'photoionization cone' of the former resembles smaller- 
scale structures in EELRs of lower luminosity Seyfert galax- 
ies; with its higher and more uniform (in terms of solid angle 
coverage) circumnuclear obscuration giving rise to shock- 
powered line emission, IRAS F23060+0505 is perhaps at a 
somewhat earlier evolutionary stage; AGN photoionization 
may become more important as the central regions are swept 
clear of obscuring material. 

Turning now to the lower luminosity objects, the high 
[OIIIJ/H/3 and [NeIII]A3869/[OII]A3727 ratios suggest that 
AGN photoionization is the dominant mechanism in IRAS 
F01217+0122. The case of IRAS F01003-2238 is more am- 
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biguous: it has a much lower [OIII]/H/3 ratio than the other 
objects and on an [OIII]/H/3 versus [NII]/Ha diagnostic di- 
agram, it falls just within the HII region area. The young 
stellar population could thus be just as important a source 
of ionization as a putative AGN. The lack of high ioniza- 
tion lines ([NeIII]A3869 is detected at just less than the 
3cr level) suggests, moreover, that the former mechanism 
may dominate. Indeed, the number of number of 05 stars 
needed to account for the extinction-corrected blue contin- 
uum could plausibly power the entire bolometric luminosity 
of the galaxy. Thus, a comparison of the roles played by 
AGN and HII region photoionization in IRAS F01217+0122 
and F01003-2238 is also consistent with their locations on 
the putative evolutionary sequence. 



8 CONCLUSIONS 

We have presented the results of an area spectroscopic study 
of the galactic-scale nebulae around four ULIRGs, focus- 
ing on the kinematics and ionization structure of the line- 
emitting gas. We accommodated our findings within the evo- 
lutionary scenario of Sanders et al. (1988a), by showing how 
the kinematics may be used to assess the merger status in a 
way which is consistent with other indications of the latter. 

The systemic velocity fields in the four objects are of 
essentially two kinds: IRAS F23060+0505 and F01003-2238 
exhibit orderly dipolar rotation, whilst IRAS F20460+1925 
and F01217+0122 show little coherent structure in the sys- 
temic together with FWHM > 1000 km s^ 1 , which we 
claim are indicative of more advanced mergers. For IRAS 
F23060+0505 and F20460+1925 we support this interpreta- 
tion with morphological and colour information derived from 
JKT and archival HST data. At lower luminosity, fits to the 
continuum spectra of IRAS F01217+0122 and F01003-2233 
are consistent with the less advanced merger having the 
younger starburst population. Furthermore, the emission- 
line spectrum of IRAS F01217+0122 clearly reveals an AGN, 
whilst IRAS F01003-2238 has an HII region spectrum and 
could be powered entirely by the starburst: we are plausi- 
bly seeing two systems at different stages along the merger- 
induced transition from starburst to AGN. 
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APPENDIX A: EMISSION-LINE FITTING 
SYSTEM ATICS 

Throughout this paper, the [OIII]+H/3 profiles from the AR- 
GUS data have been fitted with a model comprising two 
equal width gaussian velocity components, the spectral res- 
olution and signal-to-noise ratio being insufficient to jus- 
tify the fitting of anything more complex. Whilst such mod- 
els provide a good parameterisation of the blue asymmetric 
line profiles, the physical interpretation of the velocity off- 
sets and relative intensities so derived may be subject to a 
number of systematic effects. This is especially likely if, as 
found by AHM for a sample of lower- luminosity IRAS galax- 
ies, the true profile comprises a systemic component and a 
much broader, blueshifted base. 

To investigate these effects, we performed some simu- 
lations. Synthetic, noise-free, [OIII]+H/3 profiles were gen- 
erated, comprising a systemic gaussian velocity component 
(centroid Xnarrow, width a narr ow) and a second, broader, 
component with width abroad, separated from the systemic 
by Avtrue, and with a fraction rtrue of the latter's flux (the 
[OIII]/H/3 ratios of the two components were set equal). 
The profiles were then fitted with the aforementioned two- 
component, common-width model and the responses of the 
fitted parameters (A/it, &fit, Avfu and rf it ) to variations in 
the true profile parameters were investigated. The results 
are as follows: 



(1) With a narrow — 5.5A and cr broad = 8.25A (correspond- 
ing to FWHM of 500 and 1000 km s" 1 respectively, after 
correction for an instrumental resolution of 4.2 A), it was 
found that af it varies little as a function of Avtrue (remain- 
ing near 5.8A) and that A/jj follows very closely the move- 
ment of Xnarrow In contrast, the value of tfu declines from 
0.25 to 0.03 as Av tr ue varies over 900 to 100 kms" 1 and 
the behaviour of Avfu is even more unexpected: over the 
same range in Atitrue, Avfu changes from ~ 1000 to just 
~ 700 kms" 1 , as illustrated in Fig. Al. At the highest val- 
ues of Avtrue, the model underpredicts the flux between the 
[OIII]A5007 and A4959 lines. 



(2) Xfu, Avfu and r/it are insensitive to changes in 

anarrow • 

(3) As abroad is increased to 12A, it is observed that A/it 
remains practically unchanged, while r/i t and a jit decrease 
and increase by small amounts, respectively. Once again 
Avfu exhibits a strong response: as Avtrue varies over 900 
to 100 kms" 1 , Av fit ranges over ~ 1190-1050 kms" 1 . The 
flux deficit between th e [Q UI] lines becomes significantly 
more pronounced. Fig. A2 shows rju and Avfu as func- 
tions Of Avtrue- 
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Figure Al. Plots showing the variation of the fitted relative 
intensity r/it and velocity offset Avf it as functions of Av tr ue, 
for a„arrow = 5.5A, a br0 ad = 8.25A and r tr ue = 0.4 
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Figure A2. Plots showing the variation of the fitted relative 
intensity r/it and velocity offset Avf it as functions of Avtrue, 
for anarrow = 5.5A, a br oad = f2.0A and r tr ue = 0.4 



is of interest: eg. with the latter set at 12A, aj it varies from 
about 6 to 6.5 as r true varies over the above range and the 
flux deficit between the [OIII] lines remains significant. 



(4) As r t rwe is increased from 0.4 to 0.8, Xfu and Avju 
scarcely respond, and r/; t increases by roughly a factor of 2. 
The response of aju, particularly at large values of abroad, 



Fig. A3 illustrates how, even when a visually good fit 
to the complex is obtained with the common-width model, 
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Figure A3. The dots sample the synthetic [OIII]+H/3 pro- 
file with abroad = 12. OA and a nar row = 5.5A, and Av trU e = 
300kms~ . The solid line shows the best-fit two compo- 
nent, common- width model, with systemic (dashed) and 
blueshifted (dotted) components having Of it — 6.0A, sep- 
arated by Avfu — 1022 km s^ 1 . Note the great disparity 
between the model and fitted velocity differences, and be- 
tween r t rue = 0.4 and rf it = 0.085. A flux deficit between 
the [OIIIJA5007 and A4959 lines is also visible. 

the extracted parameter values can poorly represent those 
of the true profile. 
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